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Abstract 

We study localization and delocalization in a class of non-hermitean Hamil- 
tonians inspired by the problem of vortex pinning in superconductors. In various 
simplified models we are able to obtain analytic descriptions, in particular of the 
non-perturbative emergence of a forked structure (the appearance of "wings") 
in the density of states. We calculate how the localization length diverges at 
the localization-delocalization transition. We map some versions of this problem 
onto a random walker problem in two dimensions. For a certain model, we find 
an intricate structure in its density of states. 
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1 Introduction 



Non-hermitean random matrix theory 0] has been apphed recently to a number of 
interesting physical situations. An interesting issue among these physical situations is 
the study of localization-delocalization transitions in non-hermitean random Hamil- 
tonians. The earliest discussion of localization in the context of a non-hermitean 
hydrodynamics problem appears to be in 0. Later on, Hatano and Nelson have 
mapped the problem of the vortex line pinning in superconductors to a problem in- 
volving a non-hermitean random Hamiltonian. When a current is passed through a 
superconductor, vortex lines tend to drift in a direction perpendicular to the current, 
but this tendency is counteracted by impurities on which the vortex lines are pinned. 
It is expected that at some critical current the vortex lines become unpinned or delo- 
calized. In the simplest model of this problem, the physics is modelled by a quantum 
particle hopping on a ring, whose rightward (or counterclockwise) hopping amplitude 
te^ /2 is different from its leftward hopping amplitude te^^ /2. Note that ih may be 
thought of as an imaginary gauge field. On each site of the ring is a random potential 
Wi which tries to trap the particle. The Hamiltonian 

H = Ho + W (1) 

is thus the sum of the deterministic non-hermitean hopping term 

Hoij = \(e'' Si+i^j + e-"" 5i,,+i) , i,j = l,---N (2) 

(with the obvious periodic identification i + N = i oi site indices) and the hermitean 
random potential term 

Wij = Wi 6ij . (3) 

The number of sites N is understood to be tending to infinity. This problem has been 
studied by a number of authors @, ||, ^ following Hatano and Nelson. Note that 
the Hamiltonian not only breaks parity as expected but is non-hermitean, and thus 
has complex eigenvalues. It is represented by a real non-symmetric matrix, with the 
reality implying that if E is an eigenvalue, then E* is also an eigenvalue. 
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Figure 1: The spectrum of the Hamihonian in Eq. (|T]) for h = 0.5, t = 2 and N = 400 
sites. Shown here is the spectrum for one particular reahzation of site energies taken 
from a flat distribution with —2 < Wi < 2. A finite fraction of eigenvalues has clearly 
snapped onto the real axis. (The coordinate axes have been suppressed to make the 
appearance of the snapped eigenvalues clearer.) 

With no impurities {wi = 0) the Hamiltonian is immediately solvable by Bloch's 
theorem with the eigenvalues 

27rTi 

Er, = t COS ~ih), (n = 0, 1, ■ ■ ■ , iV - 1) , (4) 

tracing out an ellipse. The corresponding wave functions ipj^^ ~ exp 27rinj/N are 
obviously extended. Note that in the limit of zero non-hermiticity {h = 0) the ellipse 
collapses to a segment on the real axis as expected. With impurities present the 
spectrum can be obtained by numerical work as was done extensively by Hatano and 
Nelson and as is shown in Fig. 1. 

Two "wings" have emerged out of the two ends of the ellipse. Some eigenvalues 
have become real. Evidently, the "forks" where the two wings emerge out of the 
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ellipse represent a non-perturbative effect, and cannot be obtained by treating the 
impurities perturbatively. It is thus something of a challenge to obtain the two wings 
analytically. In this paper, we address this and other problems. 

As discussed in Section 6 of this behavior could be understood qualitatively 
by a simple example. In ordinary hermitean quantum mechanics, it is a familiar 
textbook dictum that nearby eigenvalues repel. In contrast, two nearby eigenvalues 
in the complex plane, separated along the imaginary direction, attract each other 
under a hermitean perturbation. To see this consider the 2x2 matrix 



each other, but remain on the imaginary axis. However, as soon as w > e, they snap 
onto the real axis and start repelling each other. Let us start with the ellipse in 
the absence of impurities. Near each tip of the ellipse, there is a pair of eigenvalues 
separated slightly along the imaginary direction and lying on opposite sides of the 
real axis. They attract each other and thus approach the real axis, but as soon as 
the two "friends" arrive on the real axis, they immediately repel each other (as well 
as the eigenvalue already on the real axis.) Obviously, this process repeats itself with 
the next pair of eigenvalues, and thus leads to the formation of the wings. 

Hatano and Nelson emphasized that Hq has a special property, namely that by 
a (non-unitary) gauge transformation, all the non-hermiticity can be concentrated 
on one arbitrarily chosen bond, without changing the spectrum. This leads to an 
extraordinarily simple argument that the states corresponding to complex eigenvalues 
are extended, that is, delocalized. Let Hip = Eip. Assume that ip is localized around 
some site j. We can always gauge the non-hermiticity to a link which is located 
arbitrarily far away from the site j, where |^| is exponentially small. Thus, if we 
cut the ring open at that link, the effect on the Schrodinger equation would be 
exponentially small, and would vanish completely in the limit ^ oo. It follows 
from this simple argument that if we replace the periodic boundary condition in 




with eigenvalues are 



— w'^. Thus, for w < e, the original eigenvalues ±ie attract 
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solving Hij) = Etjj by an open chain boundary condition, the locahzed part of the 
spectrum of (0) would not be affected. But for H with an open chain boundary 
condition, the gauge transformation just mentioned may be used to gauge away the 
non-hermiticity completely, meaning that the Hamiltonian is in effect hermitean with 
real eigenvalues. We thus conclude that all localized eigenstates of (|1|) correspond to 
real eigenvalues (in the large N limit.) In other words, the states corresponding to 
complex eigenvalues are extended, that is, delocalized. 

Remarkably, non-hermitean localization theory is simpler in this respect than the 
standard hermitean localization theory of Anderson and others To understand 
the localization transition, one has to study only the density of eigenvalues, or equiv- 
alently, the one-point Green's function, rather than the two-point Green's function. 
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2 The single impurity case 



In this section our philosophy is to find the simplest version of the hopping model 
described in the Introduction which we can solve exactly, but yet manages to capture 
the essential physics involved, including the non-perturbative emergence of the two 
"wings" along the real axis. The desired simplification is to replace the N random im- 
purities in by a single impurity. With this simplification we do not need to specify 
the precise form of the non-hermitean Hq beyond assuming that it is translationally 
invariant. We thus replace the Wij in by 



w 6ii6 



(5) 



where w is drawn from some probability distribution P{w). 
Our task is to calculate the averaged Green's function 



w 



) 



(6) 



of the random Hamiltonian H = Hq + W, from which we may calculate the eigenvalue 
density of H. 

Away from the location of the spectrum of Hq in the complex z plane we expand 
(P) in powers of l/{z — Hq), and thus obtain 



G{z) = Goiz) 



1 

N 



d 



1 



dz \z — Hn 



0/ 1,1 



fe=l 



1 



Ha 



0/ 1,1 



k-l 



where 



Go{z) = tr ^ 
i\ z — 



Ho 



) 



(7) 



(8) 



is the Green's function of Hq. Due to the translational invariance of Hq we have 

1 



Hn 



Gq(z) 



0/ 1,1 



and so 



Giz) = GQ{z)-j^^^py)[GQ{z)] 



k-l 



Gq{z) 



1 OGq^z) 



N dz '1-wGq(z 



w 



(9) 



Observe that the effect of the single impurity on the Green's function is of order 
as should be expected. 

We stress again that up to this point we did not adhere to any specific transla- 
tionally invariant Hq nor did we specify any particular probability distribution P{w) 
in our derivation of (^. We also note that our derivation is exact for any value of A^. 

For finite N the singularities of Go{z) are isolated simple poles located at the 
eigenvalues of Hq. The effect of the single impurity on any of these poles would be to 
move it around in an amount which depends on the typical scale r of the distribution 
P{w). Thus, ii z = zq is one of the poles of Go (-2)5 "we expect that for small values of 
the scale r, the full Green's function G{z) will also have a pole near z = zq. Thus, in 
the vicinity oi z = zq (ignoring all the other poles) we may approximate 

Substitution of (0) into (^ yields the simple result 

J - (11) 

iV Z — Zq — 

Eq. ( |lTl) is nothing but the result of first order perturbation theory (after all, we 
ignored the effect of all states other the eigenstate associated with 2:0.) 

We can now understand the non-perturbative emergence of wings under very 
general circumstances. Evidently, besides the "trivial" poles just mentioned, G{z) 
also has a pole whenever wGq{z) = 1 with w in the support of P{w). We have 

Go(^--)-/<'V-^|^ (12) 

where po{x, y) = {l/N) J2i ^{^ ~ Re Ei) 6{y — Im Ei) is the density of eigenvalues Ei 
of Hq. All we require for the following discussion is that po{x,y) = po{x, —y). This 
is true for all the Hamiltonians Hq considered in this paper. Then on the (positive) 
real axis Gq{x) is real, and decreasing outside the spectrum of Hq. (Indeed, it is 
well known that we can interpret the real and imaginary parts of Gq{z,z*) as the 
electrostatic field E = [E^^Ey) = (Re Gq, — Im Gq) generated by the charge density 
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Pq.) Thus, if the real quantity 1/w lies between Go{xedgc) and zero (where Xedge 
denotes the intersection of the edge of po with the real axis), we will have a pole on 
the real axis at some Averaging over w we thus obtain a wing on the positive 

real axis. A similar discussion can obviously be given for the negative real axis. It is 
also clear that there is no solution of wGq{z) = 1 for z outside po and away from the 
real axis. 

The probability distribution P{w) will in general depend on some set of parameters 
{vi} and for some given {rj} it is of course possible that 1/w does not lie between 
G'o(a;edge) and zero. Thus, for some critical values {r^} there will be a transition at 
which the wings, and hence the localized states associated with them, disappear. 

In the N oo limit the eigenvalues of Hq become dense and will either trace 
out a curve in the complex energy plane (analogous to the ellipse associated with 
@), or fill out a two dimensional region (as for example, in two dimensional hopping 
problems.) We now focus our attention at the hopping Hamiltonian Hq in (|^). To 
make things as simple as possible we let the parameters in (^ tend to the (maximally 
non-hermitean) limit h ^ oo and t ^ such that 



t e' 



2 



(13) 



(and obviously t e 



0). In this limit (^ changes into 



27rm 



(n = 0,l,--- 



iV-1), 



(14) 



= exp 



N 



and the ellipse associated with (^ expands into the unit circle. Furthermore, the full 
Schrodinger equation {Hq + W)ip = Eip becomes simply 



(15) 



As a consequence. 




(16) 



from which we extract the localization length of the state i/j simply as 



At the points where the wings join onto the circle the locahzation length diverges as 
expected. 

The Green's function associated with (|T^) 

1 1 

n=0 ^ TV 

may be approximated in the large N limit by 
GJz) ■ 1 





0, 


\z\ 


^) { 




\z\ 



Substituting this expression into (||) we obtain 

r 0, 1^1 < 1 

G{z) = (18) 

[ 1 + l^l>i- 

Note that independently of P{w) the disk inside the unit circle remains devoid of 
eigenvalues. To get a hold of the correction G{z) — (l/z) in the outer region it is 
instructive to carry out some explicit calculations with particular probability distri- 
butions P{w). 

As our first concrete example, we take 

p(nj) = ^[5{w - r) + 5{w + r)] (19) 

with some scale r. We then find from (flM) that 



1 / 1 1 \ 

^ ^ z 2Nz^ \z-r z + rj ^ ^ ^ ' 

If r > 1, then G{z) has two new poles on the real axis dX z = ±r, each with a residue 
1/2N. The critical value for r to induce these two new poles is = 1. The existence of 
this critical value is a non-perturbative phenomenon (though its particular numerical 
value is retrospectively not surprising at all, being set by (p!^).) From (p!7| ) we find 
that the localization length of the states associated with these poles is L{r) ~ 1/log r, 
and they thus become extended as r rc = 1. 
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Consider next the box distribution 

PH = ^e{v'-w') (21) 

for which 

The eigenvalue density p{x, y) (with z = x + iy) is related to the Green's function by 
the general relation ||^ 

p(x,y) = - Ag(^,z*) (23) 

TT OZ* 

which givesQ 

p(^, y) = ^^(y)^(^ - 1^1) ' 1^1 = 1^1 > 1 • (24) 

As in the previous example, non-perturbative effects generate a critical value for V, 
namely Vc = 1. For V > Vc, the density of eigenvalues develops two symmetric 
wings of length V — Vc such that the fraction of eigenvalues residing in these wings 
is -^(1 — ■^) . The localization length L{E) of states that reside in the wings is finite 
for \E\ > Vc and diverges logarithmically as \E\ ^ Vc- 

As our final example for this section we consider the Lorentzian distribution 

with its long tails extending to infinity. The fraction of realizations in which an 
impurity potential is stronger than the unit scale set by ( [T^ ) is (I/tt) arctan 7. We 
find that outside the unit circle 



^ ^ z Nz 



e{lm z) ^(-Im z) 



\z\ > 1 (26) 



z + 27 z — t'f 

and thus using ( p3D we find that as in the previous example, the density of eigenvalues 
develops wings along the real energy axis given by 

Pwings(x,z/) = ^ 9ix' - 1) 5{y) . (27) 



""^Here we also use the fact that in the cut plane (with the cut running along the negative real 
axis) {d/dz*)\ogz = {l/2){dx + idy)\ogz — —Tr9{—x)6{y) 
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Due to the long tails of (pSj), the wings extend to infinity, and the fraction of states 
that reside in them is (l/A^vr) arctan 7, namely, the fraction of "strong" impurity 
realizations divided by N. Note also that there is no critical value for 7, namely, the 
wings appear for all values of 7. 

Note that to obtain the analog of in the case of two impurities already involves a 
non-trivial combinatorial calculation involving Gu, Gjj, and Gij, where i and j denote 
the locations of the two impurities. (It is clear though, that the results of this section 
as they stand, are still applicable to the problem with many impurities, as long as the 
separation between impurities is larger than the localization length.) Remarkably, 
however, we show in the next two sections that taking the maximally non-hermitean 
limit ([T3|) which restrict the particle to "one way" hopping, we can readily treat the 
generic problem of many impurities. 

Finally, we stress that our derivation of (^) was not limited to one dimensional 
hopping. Given Go{z) for a translationally invariant Hq in any number of dimensions, 
G{z) as given by is still valid outside the support of the spectrum of Hq. In 
general, the spectrum of Hq in higher dimensions will fill a two dimensional region in 
the complex energy plane. In order to calculate G{z, z*) in that region, and in fact, 
to see how that region is affected by the impurity, we might have to resort to the 
method of Hermitization discussed in ^j. 
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3 Maximal non-hermiticity and many impurities: 
"One Way" models 

In the preceding section, our exact result was obtained for any value of the param- 
eter h. We may perhaps hope that we can study the many impurities problem in the 
particularly symmetric case provided by the limit (13), namely t e'* — > 2 , t —>■ 0, 



in which the ellipse of eigenvalues becomes the circle (Q). Thus, we will study the 
case of maximal non-hermiticity in which H is given by 

Hij = Si+ij + Wi Sij , i,j = 1,---N (28) 

(with the obvious periodic identification i + N = i of site indices.) We refer to this 
class of models in which the particle only hops one way, as "one way models." 
We have found a particularly simple example in which 

P{wi) = ^ [5iw, -r)+ 6iw, + r)] . (29) 

We now proceed to calculate the density of eigenvalues Ei in this "one way sign 
model" , using the master formula 

p{x,y) = (^^5(x-Rei?,) %-Imi?,)) 

i 

where we defined z = x + iy {a. careful derivation of (|30|) is given for example in 
Section 2 of |^.) It is at this point that the simplification associated with taking the 
large non-hermiticity limit may be appreciated: the determinant of z — H for H in 



(|28|) is simply 

det{z-H)=(l[{z~wA-l (31) 
\k=i ) 

(For arbitrary t and h the corresponding formula is considerably more complicated.) 
Note that ( pT]) is completely symmetric in the {w,}, and thus, for a given set of 
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site energies it is independent of the way the impurities are arranged along the ring. 
Averaging over the impurities we obtain 

(log det {z - ^)) = 2-^ f: f M log \{z - rT{z + r)^"" - ll . (32) 

In the limit N ^ oo, the binomial coefficient appearing in (|32[) is sharply peaked 
as a function of n around n ~ N/2, and thus to ffist approximation,^ 

r N N 1 

(log det - i/)) ~ log {z-r)-{z + r)- -1 . (33) 



In order to calculate the density of eigenvalues we may now insert (|33|) (and its 
counterpart for W) into (pUD, but due to the simplicity of (^) we can avoid doing 



so and simply identify the branch singularities of the right hand side of (|33|). These 
singularities clearly occur at 



2„ = ±Vr2 + e— , n = 0,l,---,A^/2, (34) 

which define the support of the density of eigenvalues. 

To summarize, we see that the original unit circle spectrum of the deterministic 
part of (^8]) is distorted by the randomness (^9]) into the curve 



= r^ + e'\ 0<e <2n (35) 
in the complex z plane. Clearly, Tc = 1 is a critical value of r. For r < 1, the curve 



(|35|) is connected, whereas for r > 1 it breaks into two disjoint symmetric lobes that 
are located to the right and to the left of the imaginary axis. The lobe on the right 
intersects the real axis at -Emin = Vr"^ — ^ and at Emax = v^rM-T. As r is decreased 
to Tc = 1 the two lobes touch at the origin and merge into a single curve when r 
becomes smaller than Tc = 1. 

In the three figures below we have plotted the curve (^5|) for three values of r in 
the different regimes r < l,r = Vc = 1 and r > 1, on top of the corresponding results 
of numerical simulations. 
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We also get the correct normalization since in this approximation: 2 ^ Nl/ ((N/2)l)'^ ~ 1. 
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-■>v ■ ■■■■ 

■ 




Figure 2: The spectrum of the Hamiltonian in Eq. (p8|) . The site energies are taken 
from p9| ) for various values of r. Figures (a), (b) and (c) correspond to r = 0.9, 1 
and 1.1, respectively. The solid curve represents Eq. (^5]). 
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The width of the scatter of the numerical resuhs around the analytical curve can 
be roughly estimated for any given Zq on (|35D by keepping (in addition to the leading 
term) all terms in the sum (Bl) with \n — {N/2)\ < Ne{zQ), where e{zo) is determined 



by steepest descent. Then the branch singularities on the right hand side of (|3^) 
would have been bounded between the curves 



(Note added: In a forthcoming publication [|T^], it is shown that the width of the 
scatter vanishes as — oo and that the spectrum of the Hamiltonian in ( PB]) is 
actually one dimensional.) 
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4 Many impurities and critical transitions 



In this section we push the analysis of the previous section further and replace (|29|) 
by a generic probability distribution P{wi). 

Using the determinant ( pT]) ^{z) = (nfcLi(-2 — "U^j)) — 1 we write the Green's 



function for the "one way" Hamiltonian (|28|) as 

Since each term in the previous equation factorizes, we finally arrive at the general 
formula 



(37) 



where in each term we average over w against the arbitrary distribution P{w). A sim- 
ple check reveals that (|37|) is consistent with (^) and (|33D . The averages {{z — w)^'^) 
may be obtained of course by deriving the generating function 

«W = (r^>. (38) 

It is worthwhile to mention here that a derivation of (^) by expanding in powers 
of the "one-way" hopping piece Hq = J2i + 1| in (EHf) has the nice feature that 



starting from any site i on the chain, the terms that contribute to the trace (even 
without taking the average) are precisely those in which the particle hopped an integer 
number of complete revolutions around the chain (which is why the averages in (|37|) 
are all raised to the power A^.) This observation may provide a physical explanation 
of why the determinant ( pT] ) is completely symmetric in the site energies wf the 
particle visits all sites equally as it hops, wherever the Wj's are. 

As a concrete application of (0) we now concentrate on the Lorentzian distribu- 
tion (p5|) P{w) = {'j/tt) (w^ + 7^)~^. We have g{z,z*) = {z + sgn (Im z))~^ and 
obviously 

((^)') = f ^. \t =^9{z,^*)f. (39) 

\z — w J yz + z7 sgn (Im z) I 



15 



Thus, using (0) we obtain the exact result 



r~'t^^*\- 9{z,z*) _ {z + i-f sgn {Im z)) 

^ ' ~ ^ — TT — 7^ ~ ~t — 7~- rf \^ — 7' ^ ' 

1 — {g{z, z*)) [z + ^7 sgn (im z)) — 1 

We would like now to calculate the density of eigenvalues. We first investigate p(x, y) 
off the real axis. Using (|3D {p{x,y) = (l/vr) {d/dz*) G{z,z*)) we find 



p{x,y) = - {z±i-f)^ ^ ^ 



TT ' " dz* {z±i-f)^ - 1 

N-l 

= mT. ^^^^ (-2 + ^7 sgn (Im z) - uJk) , (41) 

fc=0 

where Uk = e^'^*'^/^, and z ± ^7 in the first line of ( ^T]) correspond to z being in 
the upper or lower half plane, respectively. The complex eigenvalues are thus equally 
spaced along the union of two arcs of a circle of radius one. The upper arc is that part 
of a semicircular arc of a unit circle (centered at the origin) that remains in the upper 
half plane after being pushed a distance 7 downward along the imaginary axis. (The 
lower arc is of course the mirror image of the upper arc.) Each of these arcs is thus 
of length arc cos (27^ — 1) < vr and carries narc = {N/2tt) arc cos (27^ — 1) < {N/2) 
eigenvalues. This means that the arcs exist only as long as 7 < 7c = 1, which should 
be contrasted with the situation in (p7|). There, the single impurity does not perturb 
the unit circle, which persists for all values of 7. 

The rest of the eigenvalues, which do not have space to live on the arcs, must 
have snapped onto the real axis and formed "wings". The eigenvalue density along 



the real axis is generated when d/dz* hits the step-functions 6'(±Im z) in (40) and 
we find that it is given by 

1 XV w {x-i-i)^-^ 6{y) /i^ sin + sin (A^ - 1)0 

p...,(x,,) = - 6iy) = P _ cos iV0 + 1 ' 

(42) 

where we have defined a; — ^7 = p e*"^. 

In the large N limit (^) tends to a particularly simple form. It is clear that this 
form depends on whether /x < pc = 1 or /i > /ic. For /i > 1 

Pwing(x, y) = -1 ^ eip' - 1) Siy) = ^ 9{x' + f - 1) Siy) (43) 
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(which again should be contrasted with the one impurity case (0),) while for yU < 1 

Pwing 0. 

We thus see that for 7 < 7^ = 1, there are two wings that bifurcate from the arcs 
at X = ±Xc = ±v^l — 7^. Integrating over (^) we find that the fraction of eigenvalues 
that reside in the wings is (n^ings/N) = 1 — (2/7r) arc cos 7 = 1 — (l/vr) arc cos (27^—!), 
which together with the fraction of eigenvalues 2{naic/N) = (I/tt) arc cos (27^ — 1) 
that reside in the arcs, sum up exactly to 1. As 7 tends to 7c, Xc becomes smaller, 
and vanishes at 7 = 7c. At this point the two wings touch at the origin and the two 
arcs disappear completely. 

As already mentioned, there is no critical 7c in the single impurity perhaps 
might be expected. The many impurities case is dramatically different: the long tail 
of the Lorentz distribution can overwhelm the non-hermiticity for 7 > 7c, and the 
spectrum collapses to the real axis. We have carried out some numerical studies. 
For 7 away from 7c = 1, our analytic results fit the numerical data closely, but 
as 7 approaches 7c, the statistical fluctuation between different realizations of {wi} 
becomes larger and larger (at 7 = 0.9 for example, for as large as 400.) It would 
be interesting to study the character of the transition in more detail. 

It is perhaps worthwhile to remark that although in many discussions of disordered 
physics the gaussian distribution is the simplest to deal with, here it leads to a g{z) 
not given by an elementary function. 

Having derived a closed formula (Eq. (|37|)) for the Green's function of the "one- 
way" model (p^) (with pure clockwise hopping Hq = + M), we may now 
perturb it by adding a term AH = rj^i K + (with r small) which allows the 
particle to hop counter-clockwise. To first order in perturbation theory, the correction 
to the Green's function ( pT]) is now 

r 1 / 1 / 1 \' 

tr (Ho (AH 

z — w \ z — wj \ z — w J 

which reflects the fact that the particle has completed one revolution around the 
chain by performing A^ + 1 steps counter-clockwise and one step clockwise (which 
can occur anywhere along the chain.) It is thus possible to study the non-hermitean 

17 



Hamiltonian (|1]) by perturbing both from the hermitean hmit and the maximally 
non-hermitean limit. 

In 1^ we showed that by a hermitization method we can associate with every 
non-hermitean Hamiltonian H a hermitean Hamiltonian Ti. In particular, for the 
hopping Hamiltonians H discussed in this paper, the Hamiltonian Ti involved the 
hopping of a particle with a binary internal state (call it an " up" or a " down" particle) 
such that as it hops it flips its internal state. Solving the Schrodinger equation 
associated with Ti amounts to a simultaneous solution of the Schrodinger equations 
associated with the two hermitean positive Hamiltonians Hi = {z — H)[z — H)'^ and 
H2 = {z — Hy{z — H). We see that for our "one-way" models Hi and H2 involve only 
nearest-neighbor hopping, while for the non-hermitean hopping models such as ([^) or 
its generalizations described in Section (5), Hi and H2 involve next-nearest neighbor 
hopping^ 

■^This observation persists in the continuum of course. In the continuum the "one-way" Hamil- 
tonian becomes a first order differential operator and thus Hi, H2 are second order differential 
operators. The generic hopping Hamiltonian becomes a second order differential operator and thus 
Hi and H2 are fourth order differential operators. 
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5 Continuum "One Way" models 



The discrete "One Way" models of the two previous sections were solved exactly in 
the lagre N limit. We now show that their continuum counterparts are also exactly 
solvable. 

Starting with the continuum non-hermitean Schrodinger equation 

H = -{l/m)[d^ + h]'^ + W{x) (44) 

(with constant h), we reach the "One Way" limit by letting both h and m tend to 
infinity with a finite ratio h/m (which we set to 1/2.) In this limit, H in (^) turns 
into the first order non-hermitean operator 

H = -d., + Wix) , (45) 

which is the desired continuum "One Way" model. Clearly, the spectrum of (|^) may 
be solved explicitly for any W{x), once the boundary conditions are specified. For 
example, if (^5]) is defined over < x < L with periodic boundary conditions, we 
have 

1 ^ 
i'nix) = -y= exp ( / dyW{y) - E^x) 
VL { 

X 

<Pi{x) = 4f (- / dyW{y) + Enx) (46) 
VL { 

where ipn and (pn are the eigenvectors on the right and on the left, respectively. 

L 

Imposing the boundary conditions we find En = t J dyW{y) + {n being an inte- 



ger.) Since we know the spectrum {En} explicitly, we can now consider randomizing 



W{x) in ( ^31) and calculate any desired correlation function given any distribution of 
the random W{x). For example, if W{x) is drawn from the Gaussian distribution 
P[W] = (l/Z) exp [-{l/2g^) Jq dxW\x)] we readily find that the averaged density 
of eigenvalues is 

°° f 2-Kn\ I 1 \ 

p{ReE, ImE) = ( ^ 5 [ImE - -j-j 6 \ReE - - dxW{x)\ ) 
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J2 S (im^ 



27m 



(47) 



namely, each of the purely imaginary eigenvalues of Hq = —dx, is smeared along the 
real axis by the fluctuations of the real potential W{x). The spectrum is thus one 
dimensional, as in the discrete case. 
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6 Random hopping models and their mapping onto 
random walkers 



In this section we study a different class of models, in which there is no site energy, 
but the hopping is random and non-hermitean. We consider the Schrodinger equation 



with the hopping amplitudes sj and tj generated randomly according to some pre- 
scription. 

If the Hamiltonian is hermitean, then Sj = tj. This hermitean problem was studied 
some twenty years ago by Eggarter and Riedinger 0. They mapped the model onto 
a random walk problem and were able to show that all the states, except for the 
one at E = 0, are localized, and furthermore, that the localization length diverges 
like |log-E| as ^ 0. The existence of localized states is in accordance with the 
arguments of Anderson et al. f^. In contrast, the appearance of an extended state 
at precisely E = is not generic and is due to the invariance of the spectrum under 
E —E (as one can see by flipping the sign of ipi for i odd.) Thus, the extended 
state at = is unstable under any perturbation that destroys this symmetry, such 
as adding random site energy. 

Here we extend and generalize the analysis in |^ to the non-hermitean case. Our 
discussion below serves also as a review of 0, since obviously at any stage we can set 
Sj equal to tj. Dividing ( ^Sf ) by ipj and defining = tj_iipj_i/ipj we obtain 



where Rj = s*tj. This equation is of course equivalent to Schrodinger's equation (|48|) 
and allows us to solve for Aj iteratively and hence for the wave function ip. For a 
closed chain the obvious boundary condition is A^v+i = Ai. 

Performing a gauge transformation ipj —>■ Xjipj we find that we can effectively 
transform s* (Aj_|_i/Aj)s* and tj (Aj/Aj+i)tj. As perhaps expected, Rj is invari- 
ant under this transformation. We see by these considerations that the open chain 



Eipj = s*jijj+i + tj-i^Aj-i 
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and the closed chain are quite different. For an open chain, we can always use this 
gauge freedom to set all the tj equal to 1, say, with no loss of generality (which in the 
hermitean case would mean setting all the Sj equal to 1 as well, hence getting rid of 
randomness altogether.) On the other hand, for a closed chain with sites, we have 
the boundary condition ipN+i = V'l and hence the constraint Aat+i = Ai. The gauge 
transformation is not in general useful. 

We can of course invent a variety of models. Hatano and Nelson |^ considered 
s* = Tj and tj = tj e~" where tj is taken from the flat distribution P{t) = 
2^ ^^(A^ — t"^)- We have studied one particularly simple model which we call the 
"clock model" , defined by setting sj and tj equal to random phases. In some sense, this 
model is particularly attractive, in that it includes no free parameters. Numerically, 
we found that the eigenvalues in the clock model are distributed (in what appears to 
be a uniform distribution) in a disk in the complex plane, centered at the origin and 
with radius approximately equal to 7r/2. The rotational invariance of the spectrum 
results from of our freedom to multiply the random Hamiltonian by an arbitrary 
overall phase. 

Recall that the spectrum of any non-hermitean (or hermitean) hopping problem 
as defined in (^) is invariant under E —E. Furthermore, in the case of the 
clock model, the existence of an extended state at E = (provided = is in the 
spectrum, which is always the case when the number of sites is odd) survives the 
non-hermiticity. This is so because the Schrodinger equation (|48|) for E = implies 
that = = iV'i-il) and so this state is obviously extended. We expect 

the other states to be localized with an energy dependent localization length that 
diverges as E 0. These expectations are supported by the numerical studies we 
have done. 

We would like now to discuss the divergence of the localization length as E 0. 
We thus focus on (^) for small E, that is E small compared to a typical value of A^. 

In the hermitean case, we have the important observation that Rj = is real 
and positive. Also, E is real, and thus by ( ^Dj ) we can take Aj to be real (for an open 
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Figure 3: Two representative wave functions in one particular realization of the "clock 
model". The wave function in Fig. (a) corresponds to the eigenvalue closest to the 
origin E = 0.04 — O.lii. It is extended and its participation ratio is 0.75. The wave 
function in (b) corresponds to the eigenvalue farthest from origin E — —0.78 + 1.57i. 
It is localized and its participation ratio is 0.01. 



23 



chain, of course.) We see from that for E small compared to the typical scale 
of A, the quantities changes sign from site to site, and hence, as pointed out by 
Eggarter and Riedinger, it is convenient to iterate (HUf) twice and write 



( \ 


\ 1 - — 1 




1 + 



(50) 



We can now define Zj = logA,- and interpret Zj as the position of a random walker in 



the complex plane and j as time. Taking the logarithm of (pOD we obtain 

1 



Zj+2 = Zj + log + log 



E 
A, 



E{Aj-E) 



1 + 



(51) 



defining the motion of the walker. 

As noted above, for the hermitean case, Rj is real and positive and thus the 
random walker stays on the real line. In particular, Eggarter and Riedlinger took 
(log {Rj^i/Rj)) = and defined = ((log {Rj+i/Rj))'^). They made the insightful 
observation that for ~ the last term in (0) (which would otherwise be too 
difficult to treat) could be taken into account effectively as boundary conditions set 



on the walker. Consider first the —log 



1 + 



E{Aj-E) 



part of that term. When 



A • ~ 



(Rj/E) ~ a large positive number (we can always think of Aj as positive), 
the walker's position on the real line (namely, Xj = log Aj) decreases rapidly. We 
can thus effectively replace the term under consideration by a refiecting wall located 
at log(i?*/-E) (which moves off to infinity as ^ 0), where R* is a typical value of 



Rj. Consider now the remaining log 



1 - 



part of that term. This part becomes 



important when Aj = 0{E), at which point the sequence of A^ switches sign. We 
can thus think of a trap (or a "walker-eating monster") located at x ~ logE (which 
wanders off to — cxo as ^ 0.) Eggarter and Reidliger showed that the localization 
length can be related to the lifetime of the walker. In a numerical simulation we start 
with a walker being born at the wall (x ~ +oo); typically he drifts rapidly towards 
X ~ 0, where he executes a random walk in his middle age, and as soon as he drifts 
into a substantially negative x region, he rapidly approaches his death at x ~ — oo. 
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It is perhaps satisfying that in going from the hermitean problem to the non- 
hermitean problem the random walker has escaped from the one dimensional world 
and wandered off into the complex plane. (Strictly speaking, due to the properties 
of the logarithm, the walker now lives on a cylinder with circumference 2tt.) The 
equation ( |5TD describing the two dimensional walker is considerably more difficult to 
treat. One particularly simple case is given by the clock model in which the {Rj} are 
pure random phases. Thus, separating Zj = xj + iyj, we see that as long as the term 
in square brackets in (^) is close to unity, the walker wanders in the y direction, with 
its X coordinate hardly changing. 

We can estimate the localization length L{E) rather crudely by noting that in the 
(supposedly) exponential tail of the localized wave function \Aj\ ~ exp {±1 / L{E)) 
and so Xj ~ ±1/L{E). Thus we estimate 1/L{E) to be given by some average x 
coordinate of the walker. In numerical studies, we observe that indeed, with the 
walker starting at x = he drifts into a random walk around some average < x >, 
(but then eventually wanders off.) 

One message of this section is that random (non-hermitean) hopping models ap- 
pear to be considerably more complicated than random (non-hermitean) site energy 
models. Indeed, let us mention that another interesting model we have studied, which 
we call a hopping "sign model" , is a restriction of the clock model, in which sj and 
tj are randomly (and independently) equal to ±1. This restriction destroys the ro- 
tational symmetry of the spectrum of the clock model, reducing it to a four fold 
symmetry, since eigenvalues must come in quadruplets ±E, ±E* (the Hamiltonian in 
this case is real and the symmetry of the spectrum under E —E remains intact.) 
We obtained the density of eigenvalues numerically, and it exhibits a complicated 
interesting structure as shown in Fig. (4). It is an interesting challenge to calculate 
this structure analytically. 

The hopping "sign model" is the strictest restriction of the "clock model". It is 
thus interesting to investigate how the eigenvalue distribution of the "sign model" 
changes when we allow the hopping amplitudes to take on more values. For example. 
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Figure 4: The spectrum of the hopping "sign model" for a chain with 400 sites. The 
amphtudes s and t take on values ±1 with equal probability. 
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we studied numerically a Z4 model in which sj and tj took on (independently) values 
from {±l,±i} with equal probability. We found that the intricate structure in Fig. 
(4) was largely washed out towards the uniform disk distribution of the "clock model" , 
with some residual fine structure. 

At this point we mention a class of hopping models in which the s's and t's in 
(|48|) are taken from the same probability distribution P{x) such that x is always real 
and positive. Although the Hamiltonian H is non-hermitean, the eigenvalues, that 
is, the solutions of det(-E' — H) = 0, are all real for large. While surprising at first 
sight, this fact can be easily understood. The crucial observation is that Rj = sjtj 
are all real and positive. For A^ large, it does not make any difference whether the 
chain is closed or open as far as the eigenvalues are concerned, and so we can apply 
the transformation sj s'j = {Xj^i/ Xj)sj and tj t'j = (Aj/Aj+i)tj mentioned 
earlier to make s'j = t'j = ^JRj (which is achieved by choosing Aj+i/Aj = \ftjJ7j.) 
Thus the model is effectively hermitean, that is, the eigenvalues (which are gauge 
invariant, of course) are real. Note that the "gauge" transformation just mentioned 
does not preserve the localization property of the wave function ipj. For instance, in 
numerical work, it is convenient to study the participation ratio, defined byQ P{E) = 
l/{NJ2j l^il^); such that P ^ for a localized state and P ^ 1 for an extended 
state. Clearly, P is not preserved by the "gauge" transformation in general, and 
therefore, say, a localized state may be gauged into an extended state (or vice- versa), 
unless the Xj are appropriately bounded as a function of the site index j. We thus 
have to do a case by case study of how the A/s behave as a function of j to infer 
the localization properties of the original wave functions from the gauge transformed 
wave functions of the effectively hermitean problem, which are of course all localized, 
as predicted by Anderson and others. 

Note that the "sign model" mentioned in the previous paragraph and the hopping 

^In general, in addition to the eigenvector on the right corresponding to an eigenvalue E, Hip = 
Eip, there is also the eigenvector on the left, H^'^cj) = E*(j>. With the normalization condition 
'Ylij4'*j''Pj — 1; we can naturally consider the variant definition of the participation ratio given by 

p{E) = i/{NY.,\m\^)- 
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Hamiltonian Ho in both evade the defining conditions of this class of crypto-real 
models. An interesting question in mathematical physics is to calculate the fraction 
of eigenvalues escaping into the complex plane when the support of P{x) includes 
negative values of x. 

Finally, in connection with the question raised in the last paragraph, we mention 
here a class of hopping models in which the s's and fs in (^8]) are of the form sj = t+Tj 
and tj = t — Tj, where the independent random amplitudes {Tj} take values in the 
range — T < Tj <T according to some, say, even probability distribution (this model 
is a one dimensional discrete analog of the two dimensional "model I" in 0.) For 
^2 ^ rp2 g^jj clearly positive with probability one, and thus the model is "crypto- 

real". In general, the eigenvectors of the hermitean gauge transformed Hamiltonian 
would be all localized. Eigenvalues would start migrating into the compex plane only 
when < T^. 

An amusing exception to the assertion that all states of the hermitean gauge 
transformed Hamiltonian are localized for > is provided by the case in which 
the {Tj} take on values ±T with equal probabilities. In this case Rj = R — — T^ 
is actually deterministic and site independent and describes free hopping. Obviously, 
all states of the hermitean Hamiltonian are extended, regardless of the magnitude of 
disorder T. Clearly, VR, and thus all eigenvalues are real when > T^, and become 
pure imaginary for t"^ < T"^. Thus, for this particular sign distribution, the answer to 
the question raised above is that for < all eigenvalues escape into the complex 
plane (and onto the imaginary axis.) 
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7 Conclusion 



We have studied localization and delocalization in a wide class of non-hermitean 
Hamiltonians. In various simplifying limits we were able to obtain analytic expres- 
sions. Our first model @, while it involved only a single impurity, was able to 
capture the non-perturbative essence of the localization-delocalization transition. It 
was widely applicable to a generic Hq, in any spatial dimensions. For more explicit 
expressions, we had to invoke the maximally non-hermitean, or "one-way" limit, in 
which the spectrum of Hq was a circle. (It is difficult to imagine that comparably sim- 
ple expressions can be obtained when the spectrum of Hq is an ellipse.) We were then 
able to analyze the problem of infinitely many impurities, one at each site. Results 
were given for several representative probability distributions P{w) for the impurity 
potential energy w. We analyzed the phase transitions as parameters in P{w) were 
varied. 

We studied non-hermitean random hopping Hamiltonians. It has been known that 
the hermitean one dimensional random hopping problem can be mapped into a ran- 
dom walk problem on the real line. We found that in going from the hermitean to the 
non-hermitean problem, the random walker steps off the real axis and goes wandering 
off into the complex plane. We showed by numerical work that for the random phase 
or "clock model" the density of states appeared to be uniformly distributed over a 
disk, while in contrast for the "random sign" model the density of states showed a 
fascinatingly intricate structure. 

The study of non-hermitean Hamiltonians opens up a rich area for exploration. 
We can immediately think of many questions to be answered. For instance, consider 
the many body problem (see also |^.) In the "one-way" model of (^8]) the many 
(non-interacting) fermion ground state will be described by a Slater determinant 
\I' ~ detij[ipi{xj)] which in the large N limit assumes the form ni>j [e*^' — e*^^]. 
The question arises, even before we consider impurities, on what is meant by the 
lowest energy state of Ho. In other words, as we fill the system with non-interacting 
fermions, how do we order the energy levels if they are arranged in a circle? These 
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questions may perhaps be answered by coupling the system to some agent (such as 
a radiation field or a heat bath) which can exchange energy with the many body 
system. 

As another set of questions, we can ask how a non-hermitean many body sys- 
tem can be second quantized? What are some of the properties of a non-hermitean 
quantum field theory, such as a gauge theory in which the gauge potential is non- 
hermitean? In a relativistic theory, if the Hamiltonian becomes non-hermitean, does 
Poincare invariance imply that its other nine generators should become non-hermitean 
as well? 

Acknowledgements One of us (A. Zee) thanks D. Nelson for interesting him 
in this problem and for continuing discussions, P. W. Anderson and M. P. A. Fisher 
for discussions about localization, and J. Chalker for suggesting that we examine the 
Lorentz distribution. We also thank J. Miller for calling our attention to 0. Part 
of this work was done at the Institute for Advanced Study where he was the Dyson 
Distinguished Visiting Professor and he thanks the Dyson Fund for its support. This 
work was partly supported by the National Science Foundation under Grant No. 
PHY89-04035. 



30 



Appendix : Non-hermitean perturbation theory 



It is of course a trivial exercise to treat the model in (||) perturbatively in the 
strength of the random impurity. Denote the eigenvalues of Hq and H by {E^} 
and {Ef^{w)} respectively. Write Hq = S^-^ES with E the diagonal matrix with 
elements {E^} (and thus the columns of S^^ and the rows of 5* are the right and left 
eigenvectors of Hq, respectively.) Define 

WiWj E^ij + ■ ■ ■ (A.l) 

« id 
Then by simple arithmetic we obtain 

E^i = S'j^ S^i 

c— 1 c C 

U^^ - Ey 

and so on. It follows that the averaged density of eigenvalues of H is given to 0{w'^) 
by 

^ |25'(x - Re E^) 5{y - Im E^) ^ Re E^u (A.3) 
Re E^) 5{y - Im E^) ^(Re E^,f + [x ^ Re ^ Im]|| 



2 

5"(x 



where we have used only {wiWj) = w'^6ij. The result (|A.3| ) holds for all Hq (and 



for finite A^.) For the specific Hq in (^ and in the large A^ limit there are drastic 
simplifications. We have E^^i = 1/N and so Im i?^. = and J2i (Re -E^i)^ = 1/A^. 
Thus the only non-trivial quantity that comes in is 

E Ef^ii ^Nt^^ cosh {h + tkifi)) - cosh {h + ik{u)) ^"^'^^ 



where k{fi) = 2TTfi/N, fi = 0,1, ■■■ N — 1. The expression (|A.3| ) simplifies to 



p{x, 1/) = ^ E "^(^ ~ R*^ ~ E R^ Ef^ii)^iy - Im - E ^m^*) • (^-s) 
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The sum in ( [A.5| ) can be converted to an integral and evaluated in principle. For 
the "one way" model described in the text the integral is particularly simple and we 
obtain 

p{x, y) = / ^ 5 (a: - (1 + ^)cos 9^ 5 (^y - (1 - ^)sin 9^ . (A.6) 

To this order in w, the circle has been turned into an ellipse. Of course, we do not 
see any trace of the wings in perturbation theory. 
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